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ABSTRACT 

White dwarfs can be used as galactic chronometers and, therefore, provide 
important information about galactic evolution if good theoretical models of 
their cooling are available. Consequently, it is natural to wonder if all the 
sources or sinks of energy are correctly taken into account. One of these sources 
is partial differentiation of the chemical components of the white dwarf upon 
crystallization. In this paper we use a new formalism to show that if there is 
a redistribution of the elements inside the star, there is a net release of energy 
that has to be radiated away and that slows down the cooling rate of the white 
dwarf. 

Subject headings: stars: interiors — white dwarfs. 
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1. Introduction 

The evolution of white dwarfs is essentially a cooling process that lasts for ~ 10 Gyr. 
Since the study of white dwarfs allows to obtain information about the age of the Galaxy 
(Winget et al. 1987, Garcia-Berro et al. 1988a,b, Hernanz et al. 1994), it is important to 
identify all the sources of energy as well as the mechanisms that control its outflow. 

The vast majority of white dwarfs (those with masses in the range 0.45 ^ M/Mq ^ 1.1) 
are made of a mixture of carbon, oxygen and some impurities coming from the metal 
content of the parent star. The most important of these impurities is 22 Ne, which results 
from He-burning on the ashes of the CNO cycle, and reaches ~ 2% by mass in Population I 
stars. Since during the cooling process the star experiences a phase transition, it is natural 
to wonder if a change of solubility at the onset of crystallization can provide an extra source 
of energy (Schatzman 1958, Stevenson 1980, Mochkovitch 1983). 

Segretain and Chabrier (1993) and Segretain et al. (1994) computed phase diagrams 
for arbitrary binary mixtures in terms of the modern density-functional theory of freezing. 
They showed that the shape of the phase diagram was completely characterized by the 
charge ratio of the mixture, Z 1 /Z 2 . Their diagrams evolve from the spindle form for 
0.72 < Zi/Z 2 < 1, into an azeotropic form for 0.58 < Zi/Z 2 < 0.72 and finally into an 
eutectic form for Z\jZ 2 < 0.58. Using these phase diagrams, Hernanz et al. (1994) showed 
that the chemical differentiation of pure carbon-oxygen mixtures could introduce a delay of 
up to ~ 2 Gyr, depending on the adopted initial profiles of carbon and oxygen in the white 
dwarf before freezing. The settling of 22 Ne was much more spectacular since it could keep 
the white dwarf warm during several billion years, if the abundance of neon was assumed to 
be X 22 = 0.02 (Isern et al. 1991, Segretain et al. 1994). The deposition of 56 Fe, the second 
impurity in importance, turned out to be less important (Xu and Van Horn 1992, Segretain 
et al. 1994) because of its smaller abundance, X 56 = 0.001. In any case, and contrary to 
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the carbon-oxygen case, these impurities should not affect the determination of the disk age 
from white dwarfs since their abundances are expected to be negligible in the oldest objects 
(Hernanz et al. 1994). 

Although the importance of the mechanical and thermodynamical consequences of 
the solidification of alloys has been recognized in Geophysics (Loper 1978, Chen 1995), 
this is not the case in Astrophysics and very often these phenomena are either completely 
ignored or sometimes, misinterpreted. In this paper we reexamine the efficiency of 
convective mixing, which is responsible for the redistribution of the chemical elements 
during crystallization, and we provide a suitable formalism to compute the subsequent 
energy released during this process. 



Due to the spindle shape of the phase diagram of C/O mixtures, the solid formed upon 
crystallization is richer in oxygen than the liquid and therefore denser. Using the condition 
of pressure continuity, the density excess can be estimated to be: 



where P ; and P e are the ionic and electronic pressures respectively, 7 is the electron 
adiabatic index and Y e is the number of electrons per nucleon. For a 0.6 Mq white dwarf 
with equal amounts of carbon and oxygen, Sp/p ~ 1CT 4 . Therefore, as the solid core 
grows from the center of the star the lighter liquid left behind can be redistributed by 
Rayleigh- Taylor instabilities (Stevenson, 1980; Mochkovitch 1983). 

The local energy budget of the white dwarf can be written as: 



2. The physics of the crystallization process 
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where all the symbols have their usual meaning. Assume now that the white dwarf is made 
of two chemical species with atomic numbers Z and Z 1; mass numbers A and A 1 , and 
abundances by mass X and X±, respectively (X + X± — 1). In the following, the suffix 
will always refer to the heavier component. Using the First Principle of Thermodynamics, 
the right hand side of equation (2) can be written as: 

dL L = _ e _ p dV_dE 
dm dt dt 

where E is the internal energy per unit mass and V = 1/p. We can also write 

dE _ /<9£x dT /<9£x dV , dE x dX 

~dt ~ \~df)v,x ~dT + ^dV^T,x ~dt + ^dXo^T,v~dT ^ ' 

Using the elementary thermodynamic relation 

equation (2) can now be written as: 

dL dT ( dP^ dV dM M \^( dE \ Xo f(i \ 

where l s is the latent heat of crystallization and M s is the rate at which the solid core grows; 
the delta function indicates that the latent heat is released at the solidification front. Notice 
that chemical differentiation contributes to the luminosity not only through compressional 
work, which is negligible, but also through the change in the chemical abundances, which 
leads to the last term of this equation. Notice, as well, that the largest contribution to L r 
due to the change in E exactly cancels out the P dV work for any evolutionary change 
(with or without a compositional change). This is, of course, a well known result (Mestel 
1952, Kovetz and Shaviv 1975, Lamb and Van Horn 1975, D'Antona and Mazzitelli 1990) 
that can be related to the release of gravitational energy (see below) . 
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Finally, integrating (6) over the whole star, we obtain: 

f%D rfX 1 fMyyi) /OP \ dV , 

L + L u = — C v — — dm — / i — — — — dm 

Jo eft Jo ^dT)v,x dt 

dM s f M WDfd E, dX 

The first term of equation (7) is the well known contribution of the heat capacity of 
the star to the total luminosity (Mestel 1952). The second term represents the contribution 
to the luminosity due to the change of volume. It is in general small since only the thermal 
part of the electronic pressure, the ideal part of the ions and the Coulomb terms other 
than the Madelung term contribute (Kovetz and Shaviv 1975, Lamb and Van Horn 1975). 
However, when the white dwarf enters into the Debye regime, this term provides about the 
80% of the total luminosity preventing the sudden disappearence of the star (D'Antona and 
Mazzitelli 1990). The third term represents the contribution of the latent heat to the total 
luminosity at freezing. Since the latent heat of Coulomb plasmas is small (~ /cT s /nucleus, 
where T s is the solidification temperature), its contribution to the total luminosity is modest 
although not negligible. The fourth term represents the energy released by the chemical 
readjustement of the white dwarf, i.e. the release of the energy stored in the form of 
chemical potentials. This term is usually negligible in normal stars, since it is much smaller 
than the energy released by nuclear reactions, but it must be taken into account when all 
other energy sources are small. 

The last term in equation (7) can be further expanded. Consider that the crystallization 
front is at mass coordinate M s , that in a time interval St the crystallization front advances 
by SM S , and that the extension of the convective shell driven by the Rayleigh- Taylor 
instability is AM. Then the change in the chemical abundance of the mixing region (5Xg iq ) 
can be related (assuming perfect mixing, see §3 below) to the difference between the 
chemical abundances of the liquid and the solid: 

SM S (X s ° l - X l j q ) = -AM 5X l j q (8) 
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Besides, we can write 



where (dE/dX ) Ms is the partial derivative evaluated at the edge of the solid core. Taking 
into account equation (8), we can now write 



dE 



Jo W^^T ( " o) feU"\axo/ 



(10) 



where 



dE 



dX 1 AM 



M Jam ^dX^T,v dm 



;n) 



From equation (10) we can define the total energy released per gram of crystallized 
matter due to the change in chemical composition as: 



'liq\ 



(dEs i dE 
\3X~Jm s ~ 



(12) 



^dX JM s \dX > 

The square bracket is negative since (dE/dX ) is negative and essentially depends on the 
density, which monotonically decreases outwards. The internal energy per unit mass can be 
divided into the electronic and the ionic components, so we can write e g = e c + e\. For the 
sake of simplicity we are only going to use here the completely degenerate nonrelativistic 
expression for the electrons 

E e = \k iP ^Y^ (13) 

where K x = 1.004 x 10 13 (cgs units), and the ideal contribution plus the Madelung term of 
the Coulomb energy for the ions 



E^^L (1- 0.9T) 
H V2 / 



(14) 



where 3ft is the gas constant and T = T e Z 5 / 3 is the Coulomb coupling constant, with 
r e = 2.272 x 10 5 (pY c y/ 3 /T. 
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Since the contribution to the electron mole number of carbon and oxygen are very 
similar e g is dominated by the ionic contribution (see §4 below). Thus, we can obtain a more 
transparent expression for e g showing that it also corresponds to a release of gravitational 
energy. We first write 

£g = - a (xr < 15 > 



where we have introduced the parameter 



«= (16) 

\9Xo) Ms 

The ionic pressure associated to the Madelung energy is P ; = \pEi so that 

-Zaixf - X^)U^-) ^-3a^=3a 7 ^ (17) 

5Pi being the change of ionic pressure at crystallization. The last term in the above 
equation has been obtained by taking into account pressure continuity, which leads to 
SP[ = —SP C ~ — 7-P c dp/p- Finally, using the virial theorem in the form 

3(^) - - - /Wwd (18) 

P 1 A^WD %D 

where the average is taken over the white dwarf mass and where Q = — /9GM W d/-Rwd; #wd 
and -Rwd are, respectively, the white dwarf gravitational energy, surface gravity and radius 
(f3 = 6/7 for 7 = 5/3), we obtain 

e g ~ a/3 j T^-r\ 3wd^wd— ^ kg WD R WD — (19) 
\PJp) P P 

The factor k is of the order of unity except close to the surface of the white dwarf where 

p rapidly decreases. With this expression for e g , it is clearly seen that it corresponds to a 

release of gravitational energy. 
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3. The efficiency of the mixing process 

Before computing the energy released by chemical segregation it is convenient to check 
the efficiency of convective mixing. Let v crys be the propagation velocity of the solidification 
front into the C/O mixture. This velocity, which can be obtained from the models, is very 
small, ~ 0.1 cm/yr. The mass flux of carbon released by the front in the liquid phase is: 

F% ys = pv crys (X l c - X° c ) = pv CIys AX c (20) 

where X l c and Xq are the carbon mass fractions in the liquid and the solid respectively. The 
criterion for convective instability taking into account heat conduction from the convective 
eddies is (Stevenson and Salpeter 1977): 

X > ke (21) 

with 

Hp /<9Px ( \ dp 



pel ^ dT ' 



-^ s ^)Ap£) ^ 

(23) 



1 — _ (r -\\\^R. 

T dr 13 ' p dr 



k = ^ (24) 

7~cond ~t~ Tconv 

where H P is the pressure scale height, c s is the sound velocity, p is the mean molecular 
weight, (r 3 — 1) = (91ogT/(91ogp) a d, and r cond and r conv are the conductive and the 
convective characteristic times, respectively. If I is the mixing length, the convective velocity 
can be written as 

*Wv = c s (x~ kef I 2 -L (25) 

Hp 

and the characteristic times are given by: 

1 

Tconv — \ ®) 
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7~cond 



(27) 



where K T is the thermal conductivity. If the carbon mass flux released at the crystallization 
front is to be efficiently carried by convection we have 

dX c 



^cryst 



zpconv 



PV C 



dr 



I 



(28) 



together with a small superadiabaticity 



X — ke ~ 



(29) 



If the superadiabaticity is indeed small (which will have to be checked on the final results), 
the gradient of carbon mass fraction in the white dwarf is given by 



dX c 



dr 



k 



dX c 



dr 



(30) 



ad 



where 



dXr 



dr 



48 (r 3 -i)T(§) 



p,p 



(ap\ 

\9pJ Pj t 



1 dp 
p dr 



Q 



1 dp 

p dr 



ad V 

The value of Q has been computed for a 0.6 M & crystallizing C/O white dwarf and is 
typically a few 10~ 2 . Now it is possible to obtain the Peclet number, P = r cond /r c , 



(31) 



n as 



5 v crys AX C H P 



0.1 



y crys 



f AXr 



100 cm 2 s" 



p + i 3 k t q v°- lcm /y r / V 01 / V k t J V 109cm / 

where K T = 100 cm 2 s~ 1 is a typical value of the conductivity in the white dwarf. Then 
P <^ 1 and the gradient of carbon mass fraction 



(32) 



dX, 



c 



dr 



4 x 10 



-12 



^crys 



/AXr 



O.lcm/yry \ 0.1 



100 cm 2 s" 



V 



(33) 



is so small that X c varies by less than 1% in the convective region. 



We now check that the small superadiabaticity hypothesis is correct, i.e. 



X ~ ke < x 



(34) 
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We first compute \ with the simplified equation of state described in §2. For a crystallizing 
C/O white dwarf we get 



, H P ^ 3xlO«cms-+ V7 g \ 



An upper limit of the superadiabacity can be obtained by assuming that, due to the 
interaction of convection with rotation, the Rossby number 

R ^ = Vconv (36) 

UJ I 

is equal to unity, where to is the angular velocity of the white dwarf (Stevenson and Salpeter 
1976). Then, 

_ ke= (HruA 2 = 1Q _ 7 /_Hp_\ 2 / c s \ 2 (Jj_\ 
X \ c s / V10 9 cmy V3xl0 8 cms-V \10 hy 1 J 

where II is the rotation period. We therefore conclude that, even when rotation is 
considered, convection is indeed an efficient mechanism to redistribute the carbon rich fluid 
out from the crystallization front and that the liquid phase can be considered well mixed. 



4. Consequences on white dwarf cooling 

Figure 1 displays the evolution of the energy released per unit mass crystallized as the 
solidification proceeds. For the sake of simplicity we have assumed in these calculations that 
the white dwarf is made of an homogeneous 50:50 (by mass) mixture of carbon and oxygen. 

EDITOR: PLACE FIGURE |T| HERE. 

The energy released near the center of the white dwarf is e g = 3.54 x 10 13 erg/g and 
the partial contribution of electrons and ions are, respectively e = —5.00 x 10 12 erg/g and 



- 12 - 



e ; = 4.04 x 10 13 erg/g. The electron term is negative and different from zero only due to 
the mass defect of oxygen (A = 16 — 3.18 x 1CT 4 ) relative to carbon. If we had considered 
other species with a higher number of neutrons as compared with protons, as is the case of 
22 Ne or 56 Fe, the situation would have been the reverse. The total energy released during 
this process is 1.95 x 10 46 erg. 

Since, to a good approximation, the luminosity of a white dwarf can be considered to 
be a function of the temperature of its nearly isothermal core, T c , it is possible to estimate 
the delay introduced by the solidification as: 



where T c {m) is the core temperature when the crystallization front is located at m. The 
extra time to reach a luminosity log (L/L Q ) = —4.5 is 1.81 Gyr. To compute it we have used 
the same relationship between the luminosity and the core temperature as in Segretain et al. 
(1994). Of course, the total delay essentially depends on the transparency of the envelope. 
Any change in one sense or another can amplify or damp the influence of solidification in 
the cooling of white dwarfs and for the moment there are not completely reliable envelope 
models at low luminosities (Mazzitelli 1994). 



We have provided a new formulation of the thermodynamics of phase separation upon 
crystallization that proves that chemical differentiation results in a net release of heat that 
is radiated away thus delaying the cooling of the star. This extra heating is not due to 
compressional work, but to the local changes of chemical abundances. This formulation is 
very simple and is suitable for introduction in evolutionary codes. 

We have also shown that the hypothesis of perfect rehomogeneization of the liquid is 




(38) 



5. 



Conclusions 
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reasonable and can be used to simplify the problem. Notice that if this was not the case, 
this would result in a decrease of the efficiency of the redistribution process as an energy 
source but would not invalidate the result that if there is some redistribution there is a 
release of heat. 

Finally, we want to emphasize that the total delay introduced by this extra source of 
energy depends on the phase diagrams, on the initial chemical compositions and on the 
transparency of the outer envelope. Any change in these factors can produce noticeable 
changes in the final outcome. 
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Fig. 1. — Energy released by the redistribution of elements per unit of crystallized mass 
(continuous line) and the corresponding time delay in the cooling (dotted line) as a function 
of the solid mass (lower panel) and the luminosity (upper panel). 



